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Using a new numerical code with non-uniform adapted mesh, we study the changes produced in 
the global properties of neutron stars by the motion of matter in crust region during post-glitch 
epoch. Our numerical analysis shows that these changes may contribute to explain the observed 
spin-down of rotational frequency. 
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The sudden spin jumps in rotational frequencies or 
"glitches" observed in many pulsars dij, as well as the 
transient phase "post-glitch" returning to the continu- 
ous spin-down process, are the main features of the pul- 
sar evolution and provide great information sources to 
study the structure and evolution of neutron stars 0. 

The core of neutron stars is mainly composed of super- 
fluid neutrons, protons, and relativistic degenerate elec- 
trons (protons and electrons are found in a small fraction 
of the neutron abundance) . In the outer kilometer pro- 
tons are trapped in a lattice of neutron-rich nuclei (the 
"crust"). The region between nuclei is filled by neutron 
superfluid f§|. 

The crust and the core of the neutron star are strongly 
coupled and rotate with an angular velocity f2c- On the 
other hand, neutron superfluid, in order to rotate with 
an angular velocity fig , establishes an array of quantized 
vortex lines parallel to the neutron spin axis (these vortex 
lines are surrounded by a standard matter core) [3j IB| ■ It 
is thought that "glitches" represent a variable coupling 
between the core crust component and the neutron su- 
perfluid 

The continuous spin-down due to the angular momen- 
tum loss, via electromagnetic radiation caused by the 
magnetic torque, produces stresses in the crust that even- 
tually can be broken and brought to its more spherical 
equilibrium configuration. The decrease of the moment of 
inertia produces a sudden spin jump, which may explain 
small glitches [6] but not the large (giant) ones observed 
in Vela and other pulsars. 

Several models use the superfluid as a reservoir of an- 
gular momentum by the mechanism of vortex pinning 
and unpinning. Vortices can pin to the nuclei lattice in 
the inner crust Q , not allowing for a radial component in 
their velocity. As a consequence, the superfluid stores an- 
gular momentum 5] . A sudden unpinning of vortices and 
their subsequent motion in the outer direction produce 
a rapid spin-down in the superfluid. As the angular mo- 
mentum can be considered constant, this also produces a 
sudden spin-up of the core-crust component of the star, 
which gives rise to the observed glitch. 



Different models have been proposed in order to ex- 
plain the glitch phenomena: mechanic al g litches Q , ther- 
mal glitches 9], core driven glitches [iO|, glitches based 
on centrifugal buoyancy forces pd| . catastrophic unpin- 
ning [l^ , or annihilation of proton flirx tubes in the crust 
core boundary [13]. In some of these models we found a 
motion of matter away from the rotation axis. 

In some models the post-glitch transitory epoch is 
thought to correspond to a vortex repinning period. In 
order to explain the exponential and linear behavior 
found from fits of experimental data of pulsar frequency, 
Alpar et al (1984, 1993) have considered the existence 
of different regions in the crust with a different behav- 
ior in the unpinning process. They differentiate the re- 
gions with vortex motion in the glitch and those without 
vortex motion. Also, the phonon-vortex interaction has 
been considered to explain the post-glitch spin-down [l^ . 
Another model is the crust-cracking model developed by 
Franco et al (2000), with the motion of matter to higher 
latitudes in the crust induced by starquakes. Recently, 
other model for the post-glitch epoch, based on a model 
for glitches by Ruderman et al (1998), has been proposed 

In this paper, we will study the evolution of a neutron 
star in the post-glitch epoch. In order to do that we will 
obtain a family of hydrostatic equilibrium configurations 
by solving numerically the partial differential equations 
and the corresponding boundary conditions for different 
times in this epoch. We will use a mesh adapted to the 
boundary of the star to avoid the high frequency oscilla- 
tions (Gibbs phenomenon) observed near the discontinu- 
ity surfaces (in our case the star boundary) in codes with 
non-adapted meshes 18]. 

We assume that, just after the crust cracking, the neu- 
tron star reaches a quasi-equilibrium configuration with a 
shape corresponding to the angular velocity at the glitch 
epoch (the adapting process is almost instantaneous, as 
experimental observations indicate). From this moment 
and during the post-glitch (which may take from weeks to 
months, the time needed to achieve a situation allowing 
for the formation of a new solid crust), the neutron star 
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FIG. 1: Frequency evolution in a typical large glitch in a 
pulsar of the Vela type. 



takes its quasi-equilibrium shape, in almost an instan- 
taneous manner, following the evolution of the rotation 
frequency and without the constraining stresses of the 
crust, which now is broken. 

Note that in the problem there are two different chara- 
teristic times (the rotation period and the duration of the 
post-glitch epoch). 

Our model of neutron star is composed by two regions; 
The core and the crust region. The star boundary (a sur- 
face of constant gravitational plus centrifugal potential) 
is completely free (in general, it does not correspond to 
an ellipsoid) and it changes with time. We will make the 
following approximations: the core can be described by 
a constant density perfect fluid (given the results for the 
density profile obtained in the studies that use physical 
equations of state this seems to be a good first approx- 
imation), the crust region can be described by a sur- 
face density shell corresponding to a jump in the nor- 
mal derivative of the gravitational potential on the star's 
boundary (the crust is a thin layer in the outer part of 
the neutron star), and the angular velocities of the core 
and the the crust are equal. Also, we impose axial and 
equatorial symmetry in our Newtonian study (a general 
relativistic version is under development). 

We assume that in the post-glitch epoch there is a 
motion of matter in the crust region; in our model this 
corresponds to changes in the matter density in the shell 
found in the boundary of the star which represents this 
re gion . The stratification expected in the crust region 
[sl Il7| supports this hypothesis. This matter motion may 
represent a matter motion related with vortex motion in 
the repinning process, which is favored near the poles of 
the star Q , or crust plates motion to the magnetic poles 
p^ . or magnetized patches motion Q. In our study we 
determine a matter motion in the crust region compati- 
ble, inside our framework, with the frequency evolution 
observed (Fig. [T|). This motion of matter in the crust 
region produces changes in the globlal properties of the 



neutron star. We will analyze this effect and its influence 
on glitches description. 

For each time, in the post-glitch epoch, we will obtain 
a matter distribution in the crust region as well as the 
corresponding neutron star hydrostatic equilibrium con- 
figuration. The time evolution of the matter distribution 
in the crust region and the rest of the properties of the 
neutron star, are determined by fitting the rotation fre- 
quency and the angular momentum of the calculated con- 
figurations to the values observed and predicted, respec- 
tively, for these quantities during the post-glitch epoch. 

Any hydrostatic equilibrium configuration of the fam- 
ily is obtained by solving the Euler and Poisson equa- 
tions in the interior of the neutron star and the Laplace 
equation in the exterior region. The star boundary is 
determined by the implicit equation p(r, 0) = 0, where 
p(r, 9) is the pressure; note that the boundary has to 
be found at the same time that the partial differential 
equations are solved (i.e., it is a free boundary prob- 
lem). We have to impose the following boundary con- 
ditions on the gravitational potential U; Uin\s = Uout\s 
and dn{Uin — Uout)\s = o'(^), where a{9) describes the 
surface density profile in the crust region and 9„ means 
normal derivative to the surface S. 

To solve this problem we have developed a numeri- 
cal code based on the program CADSOL [l^l (a program 
used for the numerical solution of elliptic partial differen- 
tial equations by Newton-Raphson methods, for instance, 
it has been successfully used to study Einstein- Yang-Mills 
fields [10] ). Our code uses an iterative process to cal- 
culate the surface of the star and the solutions of the 
equations on a boundary adapted non-uniform mesh. In 
each iteration the star boundary is moved and the mesh 
re-adapted to the new boundary; the Poisson and Laplace 
equations are solved with the new boundary conditions 
and then the Euler equation determines a new boundary. 
The process is repeated up to the moment that all the 
equations are verified with a given precision (lO^'^). 

Now, let us consider the spin frequency time depen- 
dence in a typical glitch. Fig. 1. The post-glitch epoch 
is fitted using an exponential dependence on top of a 
linear dependence [1] (sometimes better fits are obtained 
with several exponential terms; here we use a model with 
just one exponential term for simplicity). Then the post- 
glitch frequency can be described by the following ex- 
pression 

Kt) = J^o(t) + Ai.o[Qe-*/" + l-Q], (1) 

where i'o(t) = at + b. The meaning of the parameters 
appears in Fig. 1. The actual values of the parameters 
are fixed by the best fit of the experimental timing data. 
Here we use a simulation with the following parameters 

Q = 0.5 , T = 100 days , Aj/q = 3.35894519466 10^^, 
a = -1.34623296 10"*^ , b = 11.1964839822. (2) 

These parameters are of the order of magnitude of those 
obtained in a real large glitch. Note that the spin jump 



3 





20 40 60 80 100 
t(d) 



FIG. 2: Crust-superfluid density parameter and angular mo- 
mentum transfer evolution. 



FIG. 3: Oblaticity and equatorial radius evolution. 



poles (7r/12 in our case), with relative weights 1 : 0.01, 



value is slightly over a typical one for a large glitch. This 
is due to the fact that even when the precision obtained 
in our numerical calculations is excellent, we are in the 
threshold to fit the exponential data of a real glitch. 
Some changes in our routines are under development to 
improve the precision and avoid this limitation. 

We concentrate our attention on the one hundred days 
just after the glitch. First, we have to calculate the con- 
figuration at the glitch epoch {t ~ 0) after the spin jump. 
This configuration is characterized by the spin frequency, 
core density, crust-superfluid density (for i = 0, in a con- 
servative manner, we use uniform density), and core and 
crust-superfluid mass. The shape, angular momentum, 
and the rest of the physical properties of the neutron 
star are then obtained. 

The continuous loss of angular momentum, due to the 
electromagnetic torque, has also to be considered in the 
post-glitch epoch. Here, we assume, in a simplified form, 
that this angular momentum loss is proportional to the 
linear part of the spin frequency fit (implicitly, we are as- 
suming that in average the linear changes in the angular 
velocity of the neutron star is due, in first approximation, 
to the electromagnetic radiation; in general, other effects 
can contribute to this linear part [5|). The constant of 
proportionality is obtained from the angular momentum 
we calculate for the configuration at the glitch epoch af- 
ter the spin jump. 

Now, in the post-glitch epoch, we have the spin fre- 
quency and the estimated angular momentum of the neu- 
tron star for each time t. In order to generate our model 
we must indicate the form of the crust-superfluid region 
surface mass density. The effect we want to show is 
not very sensitive to the particular choice of the den- 
sity profile (our model allows a completely general form 
for the density profile). We choose a profile for the crust- 
superfiuid density consisting of a uniform term plus a 
Gaussian distribution centered around a latitude near the 
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Here K and 6 do not depend on 9. (Remark: Concrete 
values 0.01 and 7r/12 have no special influence on the re- 
sults described below). We also assume constant mass 
of the neutron star, constant core density, and constant 
ratio between the crust-superfluid mass Ms and the core 
mass Mc, ^ = 0.06. Then the observed spin frequency 
and the predicted angular momentum of the configura- 
tion determine the shape of the star and the actual crust- 
superfiuid region density profile. 

In order to visualize the physical properties for the 
configurations in the family we use two parameters S an 
o = -^/l — (i?p/i?e)^ (oblaticity). Note that this param- 
eter o will not caracterize completely the shape of the 
star, because it is not spheroidal. 

To obtain the values of these two parameters for each 
time t we calculate curves of constant angular momen- 
tum in the two-dimensional parameter space o — 6. In 
order to do that, we choose a value for S and then we 
move the oblaticity to compute an hydrostatic equilib- 
rium configuration with the angular momentum of the 
neutron star at the time t. The calculated angular ve- 
locity for this configuration will not be in general the 
actual one of the neutron star at this time. Owing to 
that, we have to change the value of 6, and repeating 
the previous algorithm obtain another configuration. All 
the process continues by interpolation method up to the 
moment the calculated frequency and angular momen- 
tum coincide with those of the neutron star at this time 
t, within a relative error below 10^^. We must repeat 
this iterative scheme for several values of t during this 
post-glitch epoch. 

The results for a neutron star with total mass M = 
IAMq and initial equatorial radius i?o = 12 km are pre- 
sented in Figs. 2 and 3. In Fig. 2.b we plot the evolution 



4 



of the reciprocal of the density parameter S. We observe 
that the pulsar spin-down, in the post-glitch epoch, rep- 
resents an increase of 1/S, which means an increase of 
the crust-superfluid density near the poles (matter ap- 
proaching the axis). In principle, this sounds contra- 
dictory because a motion of matter to the axis should 
decrease the moment of inertia and then a spin-up is ex- 
pected. The paradox is solved if we consider the reaction 
of the neutron star to this motion; if we look at the re- 
sults presented in Figs 3. a and 3.b, we observe that the 
star reacts increasing its oblaticity and equatorial radius, 
i.e., moving mass away from the axis. This fact compen- 
sates the motion of matter to the axis in the crust region, 
which gives rise to the final result of an increase of the 
total moment of inertia and the spin-down of the pul- 
sar. The transfer of angular momentum from the crust 
region Jc to the core is presented in Fig. 2. a. The effect 
on the frequency of this transfer is overbalanced by the 
increased moment of inertia. 

We have now a new and unexpected effect; the mo- 
tion of matter to the axis in the crust region (due to 
crust plates motion, vortex repinning or any other ori- 
gin) produces a reduction in the pulsar frequency due to 
the neutron star reaction. In general, this effect will col- 
laborate with the rest of effects proposed in the literature 
to give the final spin-down observed. 

Also, we would like to mention that the same effect has 
been observed in a constant angular momentum family of 
configurations, assuming a motion of matter to the equa- 



tor in the crust-superfluid shell. Applying this result to 
the spin jump in the glitch epoch, which can be con- 
sidered as a constant angular momentum phenomenon, 
we observe that the neutron star reacts decreasing its 
oblaticity and approaching mass to the axis. Then the 
total moment of inertia is reduced and the result is a spin 
up of the neutron star. 

Our model have some limitations such as, constant 
density for the core and the Newtonian framework. How- 
ever, we presents for the first time a global study of neu- 
tron star behaviour in the post-glitch epoch which takes 
into acc:ouiit the crust region. New numerical codes are 
under development to improve the model and we are very 
confident that this new effect will remain in the scene, col- 
laborating with the rest of known phenomena to describe 
glitches, allowing us to have a more complete picture of 
pulsar's evolution. 

We wish to thank J. Kunz and B. Kleihaus for getting 
us started on CADSOL. 
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